Possible Jahn- Teller effect in Si-inverse layers 
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Jahn- Teller effect in bivalley Si(100) MOSFET under conditions of quantum Hall effect at integer 
filling factors v — 1,2,3 is studied. This system is described by SU(4) hidden symmetry. At 
v — 2 static and dynamic lattice deformation creates an easy-plane anisotropy and antiferromagnetic 
exchange and lifts the valley degeneracy. At v — 1,3 Coulomb interaction is essential to produce 
weak easy-plane anisotropy. Three phases: ferromagnetic, canted antiferromagnetic and spin-singlet, 
have been found. Anisotropy energy of charged skyrmion excitation in every phase is found. 

PACS numbers: 71.70.Ej, 73.43.-f, 73.43.Nq 
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In a symmetric configuration of atoms the electronic 
state could be degenerate due to a high symmetry of the 
Hamiltonian. H. A. Jahn and E.Teller proved that such 
a symmetric configuration is unstable to a deformation 
lowering the overall symmetry of the Hamiltonian. This 
Jahn- Teller effect is well accounted for in molecules and 
crystals! In 2D electron gas (2DEG) under the condi- 
tions of quantum Hall effect the kinetic energy of elec- 
trons is quenched. Electron states become macroscopi- 
cally degenerate forming Landau levels. Interaction may 
remove some of this degeneracy. Nevertheless at integer 
fillings one finds the global spin- valley degeneracy. In this 
work we raise a question whether a lattice deformation - 
Jahn- Teller effect - lifts this global valley degeneracy at 
integer fillings. 

Theory of 2DEG at integer v usesi Hartree-Fock ap- 
proximation in the limit of small parameter: E c /Tioj Ci 
where E c = e 2 /Ijj is the energy of the Coulomb interac- 
tion and ui c is the frequency of the cyclotron resonance. 
At v — 1 theory predicts a ferromagnetic ground state 
with degenerate uniform spin orientation. The elemen- 
tary excitations are electron-hole pairs or neutral excL 
tons, which correspond to gapless and noninteractingtl 
spin- waves for vanishing momentum and vanishing effec- 
tive g-factor. In the limit of the large exciton momentum 
the electron and the hole become independent charged 
excitations. A special topological spin texture in 2D fer- 
romagnet called a skyrrnuonu has a unit charge and half 
the quasiparticle energyu. In this paper we study such 
textures in the presence of Jahn- Teller effect. 

In the bulk Si there is a six-fold valley degeneracy, 
which is reduced to two-fold for (100) orientation of 
2D-plane du&-.to a quantization of electron motion in 
quantum wellD. Even for orientation (111) of 2D-plane 
with six equivalent valleys only two valleys connected 
by the time reversal symmetry are occupied by electrons 
in strong magnetic fieldsQ. In this paper we study the 
effects related to this valley degeneracy and show that 
such bivalley system has close analogy to the bilayer sys- 
tems. The bilayer setup has been extensively studied 
theoretically at v = luB and at v = 2.E3 The layer de- 
generacy of electron states can be described by a four 
component spinor in combined spin-layer space. In the 
worksOy the total 2DEG Hamiltonian was subdivided 



into a symmetric part and a small anisotropic part which 
reduces the symmetry. Symmetric part is invariant un- 
der a hidden SU(4) group of transformations in the com- 
bined spin and valley Fock space which augment apparent 
group SU(2)<g>SU(2) of separate transformations in spin 
and valley spaces. 

Electron annihilation operator can be expanded using 
one electron orbital functions: 



(i) 



where a is the spin index, z is the perpendicular to 
2D plane coordinate and r is in-plane coordinate vec- 
tor. The index r = ±1 numerates two valley. Valley 
wave functions of an electron in quantum well of Si in- 
verse layer: x( z ) exp(±iQz/2), are normalized and al- 
most orthogonal for smooth real xi z ) with the negligible 
overlap J x 2 ( z ) exp(iQz)dz. Q is the shortest distance 
between the valley minima in reciprocal space, it's ap- 
proximately equal to 2/asu with asi being the lattice 
constant, Vw(fO is the wave function of in-plane mo- 
tion of an electron in valley t and it constitutes a four 
component spinor. Electrons in the system will strongly 
interact with phonons with momentum ±Q, giving rise 
to scattering processes from one valley to another. The 
Jahn- Teller effect (JTE) in Si-inverse layers corresponds 
to displacement of silicium atoms in z-direction. The new 
equilibrium is determined from the balance of electron- 
phonon and elastic energies. Strong magnetic field is es- 
sential for JTE. In this paper we study this phenomenon 
in details. 



Hamiltonian of electron gas on Si interface 

The Hamiltonian of Si inverse layer describes electron 
system interacting with phonons: 



H = H P + H„- 
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For narrow quantum well, 2DEG Hamiltonian in mag- 
netic field can be expressed in terms of spinor ^> QT (f): 
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where g is the conduction band gyromagnet ratio, \xb is 
the Bohr's magneton, B is the uniform magnetic field in 
z-direction. We neglect small terms with different val- 
ley indices in the electron density operator. We use the 
simplest Debay model for phonons with large momentum 
~ ±Q 



H ph [u z \ = — 



{d t u z f + c 2 
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d 3 p, (4) 



where u z (p) is the lattice displacement polarized in z- 
direction. We consider only the valley-mixing electron- 
phonon interaction: 



He-ph = 6^ rPi T ,(r)T± T2 ip an (r) x 
r X 2 (z)e ± ^ z u z (p)dz] d 2 r, 



(5) 



with O, k and psi being the deformation potential, the 
dielectric constant and the density of Si. r M and a^, 
with fj, — x, y and z, are the Pauli matrices in valley 
and spin space and t ± = (t x ± ir v )/2. We choose the 
units where: h = 1, e = c, B = 1 and the distance is 
measured in magnetic length Ih = 1. ^4(/?) = (0, Bx) is 
the vector-potential in Landau gauge. 

In quantizing perpendicular magnetic field we expand 
the electron operator over Landau orbital states 4> n , P {r)'- 



np,ar ? 



(6) 
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where ar and c„ PiQr are electron creation and anni- 
hilation operators in r valley with spin a. n numerates 
Landau levels and continuous parameter p in the Landau 
gauge specifies states inside one Landau level. 

We assume the quantum Hall effect conditions where 
electrons are confined to the lowest Landau level n = 
in the ground state. The electronic part of the Hamilto- 
nian contains the Coulomb and the kinetic energy on the 
lowest Landau level: 

^ = 2^E4ar^r + \ J {q)N {q)N {-q) , 

(7) 

where V(q) = 2ire 2 / nq is 2D Coulomb interaction and 
the electron density operator is: 

N(q) = c UTCp-qy <*t exp [-iq x (jp - q y /2) - q 2 /A] . 
p 

Unitary transformations of electron operators c aTltP — 
U aTl ,0T2 c pPT 2 m the combined spin and valley space leave 
the Hamiltonian (Q) invariant, where U is a matrix from 
SU(4) Lee group. For Landau level filling factor v = 



1,2,3 we assume that the ground state is uniform over 
p-orbitals with electrons of spin a% and valley Tj fill every 
orbital p of the lowest Landau Level: 
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One can easily check by inspection that any such wave- 
function represents an eigen function of the H sym 
([?]) . The state (|^) is degenerate and a set of related eigen 
states can be generated by applying uniform rotations U. 

The remaining terms in the Hamiltonian (^) consist 
from a valley splitting term due to singularity of 
the well potential on Si/Si02 interfaces, Zeeman term, 
electron-phonon interaction and the phonon energy: 



H" 
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where t is the phenomenological valley-splitting constant. 
This part breaks down the SU(4) symmetry and lifts de- 
generacy of the eigen states of SU(4)-symmetric Hamil- 
tonian (g). The splitting of energy levels is determined 
by matrix elements of weak anisotropy Hamiltonian ( |l0| ) 
projected onto a linear space of the symmetric Hamilto- 
nian (0) level degeneracy. There is no renormalizations 
of the anisotropy Hamiltonian parameters QIC ) due to 
electron-electron interaction in symmetric part (]]). Thus 
the Hartree-Fock approach is valid up to 0(H an /u>h) 
terms for the v — 1,2,3 cases. 

The total Hamiltonian (0,^O|) can be treated as in the 
theory of magnetism. There exists a local order parame- 
ter matrix Q(r) and the total bivalley energy can be ex- 
pressed in terms of this order parameter. The energy of 
isotropic part (|?]) will be expanded in powers of spatial 
variations of order parameter VQ(r) due to Goldstonc 
theorem, whereas in the anisotropy parts we retain only 
homogeneous order parameter Q. 



SU(4) Symmetric Case 

Non-homogeneous state can be generated from the 
ground state ([|) by slow rotations, with the effective ac- 
tion being dependent on rotation matrix U(t,r): 

S[U] = -i Tr log J VipV^ exp (i J £d?j (11) 

where the spin- valley symmetric Lagrangian of 2DEG is: 
d 



C 



a/3 



- J V£ ««,,J^ d 2 f+ \ J2 V(q)N(q)N(-q), (12) 



3 



H 



H 



a 




FIG. 1: Second order a) and the anomalous first order b) 
diagrams. Solid and wavy line represent the electron Green 
functions and Coulomb interaction. 



where fi* = -i U^d t U, O = -i WVU . We treat non- 
homogeneous matrix U as a classical filed. Therefore 
the effective action ( |i"l| ) describes a macroscopic motion 
of the corresponding spin-texture. The essential points 
here is the slow rotation of U on the microscopic scale 
given by magnetic length and the use of gradient expan- 
sion assuming locally ferromagnetic ground state. The 
effective action (^Tj) for multivalley was found in Ref.tl 
We derive it following the method of Ref.O. 

We consider three filling factors of uniform ground 
state (§): v = 1 with +1), v = 2 with *(f +1, J. -1) 
and v = 3 reduces to the case v = 1 under the electron- 
hole transformation. Reference state is alternatively 
given by an electron occupation matrix: N aTl ^ T2 — 
(1,0,0,0), in the case v = 1, and N aTlt p T2 — (1,1,0,0), 
in the case v — 2. Electron Green function on the lowest 
Landau level for the reference ground state of Hamilto- 
nian rtl2n with Q = reads: 



Ca ri ,/5r 2 (e,p) 



N, 



| (1 - N) aTli [3 T 

e + E Q + u-iO e + u + iO 



UT\ ,/3r 2 
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where u = —Eq/2 is the chemical potential and exchange 
constants in the limit of vanishing thickness of electron 
layer are: 



Eq — 2E\ — 
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The effective action is 5[f2] — So[£l) + S2 where 
S [Q] = iTrlog(G/G ) whereas S^fi] is represented by 
two diagrams on Fig.jjJ The first order perturbation cor- 
rection to the Green function: SG — G — Go, contains aa. 
electron propagation on the first excited Landau Leveltfl 
with the Green function: 



C'aT 1 ,Pr 2 ( e ) ~ 



Nan ,0T 



+ 



(1-N) 



€ — LU C + E\ + ^ + ^0 6 — UJ C -\- fJL -\- i0 

(15) 

The two terms of the Hamiltonian depend on the gra- 
dient matrix field: 



«/3tiT2 



(16) 



ipf3 T2 d 2 r, (17) 



where 0± = fl y i£l x (0_ = and the differential 
operators IT± shift an electron between the adjacent Lan- 



dau Levels: n_4> np (r) = \/2n(j} n ^i p (r), U + <fi 

n—lp (r) 



2n<j} np {r), though only the n = and n — 1 Landau 



Level states are relevant for our problem. An expansion 
of the 2DEG action up to the second power of the Hamil- 
tonian (|l6|) gives: 

S„[fi] - *Tr (ITiGo) + |Tt {H X G H X G ) + iTr (H 2 G Q ) . 

(18) 

The Haxiree-Fock diagram on Figjlja has been calculated 
in Refo whereas the anomalous diagram on Figjljb is 
calculated in the Appendix A. The resulting effective La- 
grangian is: 

- ^Tr[Nn + (t,r)(i- N)^(t,r)])^-, (19) 

where = -i Tr (NU^d^U). The matrices N and i-N 
are projector operators onto the physical rotations that 
form a sub-set of SU(4) Lee group. The matrix field 
il M can be expanded in the basis of fifteen generators of 
SU(4) group: {F'}, with I = 1..15. They form two com- 
plementary sets: the first even set includes those gener- 
ators that do commute with N, whereas the second odd 
set includes the remaining generators. Generators of the 
even set constitute an algebra itself. This algebra has 
a normal Abelian subalgebra formed by a single trace- 
less generator: N — v/A. A Lee group built around the 
even set of generators is called a stabilizer sub-group S 
of SU(4) group. The odd set always contains an even 
number of generators: eight in the case of v = 2 and six 
in the case v = 1, 3. 

The Hamiltonian is invariant under the time reversal 
symmetry. The time reversal operation transforms ro- 
tation matrix as U ~ * U* , gradient field as ft — > — f2 T 
and also inverts the magnetic field B z . Accordingly we 
rewrite the energy in Lagrangian ( [l9| ) in time reversal 
symmetric way: 



E eff [n] 



~2^ 



£1 



Tr (NSl^l - N)n, 



E ( 



- -y sgn(B z )e MI ,V ^ 



(20) 



where we used identity e M „Tr (fl^fl^N) — ie^V M f2^. It 
follows immediately that E e ff > 0. The first term is the 
gradient energy whereas the second term is proportional 
to the topological index of the spin texture: 



Q = 



d 2 f 

X7 O 2 _ - 7 



(21) 
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where Z is integer. The states with 2 ^ are called 
skyrmions in the name of T.H.R.Skyrme who first con- 
sidered such textures. The case Q = ±1 corresponds to 
the simplest spin skyrmions in the first valley which can 
be rotated by a SU(4) matrix to become a general bivalley 
skyrmion. The spin stiffness in E e f frffiOl) coincide identi- 
cally with that of the one- valley caseEJ/which means that 
the bivalley skyrmion energy is the same as found for one 
valley. Also the charge in the bivalley skyrmion core is 
distributed over the two valleys with long divergent tails: 



n±i(f)~±l/ (R 2 



if one neglects anisotropy. But 



the total charge of the two valleys follows a convergent 
distribution identical with a charge density in the one 
valley case 



n(r) 



2tt 



R 2 



TT{R 2 +r 2 Y 



(22) 



where R is the radius of the skyrmion core. 

Expectation value of any operator A in the ground 
state (A) — Tt(AQ) can be expressed in terms of the 
order parameter matrix: 



Q{r) = U(r)NU + (r). 



(23) 



Obviously rotations from the little sub-group S leave the 
order parameter intact. Thus, rotations in ( p3f ) can be 
restricted to a physical space of the bivalley 2DEG - com- 
plex Grassmannian manifold: U(4)/U{v)®U(A—v), We 
rewrite (EGl) in terms of Q: 



Tl (V^VO)^-^Q. (24) 



Eeff[Q] = ^ 



In this representation the topological index Q appears as 
an index of a map of 2D plane on the coset space of the 
order parameter. The rule (|2l| ) is a consequence of a well 
known homotopy group identity: 



7T 2 



C/(4) 



U(v)®U(4-v) 



Z. 



(25) 



This identity is proven in Appendix B. 

Varying the effective Lagrangian consisting of the first 
kinetic term in ( |l9| ) (rewritten as the Wess-Zumino 
termEj) and the energy ( |24| ) we find the matrix Landau- 
Lifshitz equation: 



Q dtQ 



= T v Q - 



It describes three (y = 1) and four (v 
'spin-valley'-waves with dispersion u>(q) 



(26) 

2) degenerate 
E ig 2 /2. 



local electron density. We illustrate construction of elec- 
tron density in the case v = 1 with fixed spin of electron. 
Electron annihilation operator (|l|) is expanded in basis 
of two mutually orthogonal valley wave function r = ±1. 
The density operator h(p) — ^ {p)%jj{p) has a diagonal 
element uniform over 2D plane: ^(z)(^(r)Vv(r)) = 
vx 2 {z)J2-k1 2 h ^ where average is taken over the ground 



state (|9|), and the term oscillating in z direction: Sh(p) — 
X 2 (z)(y>ji (r)ip-i(r)e lQz + h.c). We use electronic order 
parameter (p3) to express the average oscillating density 
as: 



6n(p) = 



2nl 2 H 



Tr(Q 



T+)e iQz 



Tt{Qt- 



-iQz 



(27) 



Diagonal density term couples to phonons with in-plane 
momentum q± = and can be neglected in thermody- 
namic limit. Oscillatory part of the average electron den- 
sity allows for lowering the energy by static lattice de- 
formation also oscillating with wavevector Q inside the 
quantum well: 



e 



:Sn(p). 



(28) 



This deformation is a phonon 'condensate' or the Jahn 
Teller effect. The corresponding energy gain is found 
from the minimization of H e - p h + H p h as: 



-NGTy{Qt + )Ty(Qt~). 



(29) 



where J\f — AeB /2nhc is the number degeneracy of the 
Landau level and G = O 2 J x 4 (z)dz /2iTl 2 H psiC 2 is the 
strength of electron-phonon interaction. Integral is the 
inverse width of the 2D layer in z-direction. 

The static deformation gives new equilibrium position 
of the lattice. Assuming static deformation small we ne- 
glect the change of phonon spectrum due to anharmonic 
effects. But we take into account the dynamical phonons 
which produce a polaronic effect in the second order of 
electron-phonon interaction. Expanding the phonon field 
around new equilibrium: u z (p) = u z st {p) + Su z (p), we ob- 
tain: 



He-ph = O / V4Ti(^r 2 Vw 2 (r) x 

X 2 (z)e ± ^ z ^ z 6u z (p) dz)d 2 r, (30) 



where the lattice deformation 



2Vp Sl u>(q) 



(b^ + bie-^ 



Jahn- Teller effect 

The Jahn- Teller effect is the deformation of a lattice 
lowering H e ^ p h + H v h energy. Phonons interact with the 



can be expanded in phonon creation and annihilation op- 
erators frt and bq. In order to find the ground state energy 
we sum up diagrams for the thermodynamic potential 
expanded in powers of (assumed to be) weak electron- 
phonon and Coulomb interactions using the Matsubara 



5 




FIG. 2: Two diagrams that gives easy-plane anisotropy in 
the case v = 1. Solid, dashed and wavy lines are the electron, 
phonon propagation and the Coulomb interaction 



methodt3. We express the energy in terms of the order 
parameter matrix Q. 

Electron's Green function on the lowest Landau level 
can be expressed in terms of the order parameter Q: 



G°(e n , P ) = 



Q 



l-Q 



E /2 



l£ ri 



E /2 



(31) 



We attribute the coordinate dependent electron wave- 
functions to the interaction vertices. The phonon's 
Green's function ist3: 



D{w n ,q) 



(32) 



At small transfered momentum the Coulomb line is 
much larger then the phonon propagator. In the op- 
posite case of large transfered momentum with ±Q z- 
component we neglect the Coulomb interaction. There- 
fore in the Coulomb vertex the valley index is conserved 
as well as the spin index. Due to condition Q(l — Q) = 
all one loop diagrams with Coulomb lines only vanish. 

The expression for the electron-phonon vertex in Lan- 
dau gauge reads: 



g{pp',q) 



e 



V. 



PSiC 



u P,Qy+P 



,-9i/4-9x(p+p')/2 A/ 2 



x iz - Q, 



where q is the phonon momentum and the envelope func- 
tion x 2 {lz) = J x 2 (z)e* 9zZ dz has characteristic width I -1 
in momentum space. The polaron contribution to the en- 
ergy assuming weak Coulomb interaction: Eq/lu(Q) — > 0, 
is: 



^e-ph 



Tr(G (e„,p)T + G (£„+^„,p')^) 

e n uj n pp' 

Y J 9 2 {PP',o)D{uj n ,q)=AfG(TT{QT+QT-) - ^33) 



electron energy and there is no anisotropy operator for 
a single electron because of the Pauli matrix identity: 
rf = 1. 

Coulomb interaction creates the easy-plane anisotropy 
in the case v = 1 and the essential diagrams are shown 
on Fig.(||). The dashed line represents the 2D Coulomb 
potential V(q) whereas the solid and wavy lines repre- 
sent the electron (31) and phonon (|3|) propagators. Di- 
rect calculation neglecting the dependence of uj(Q) on 
in-plane momentum q± shows that the sum of two dia- 
gram on Fig. (j^) is 



E, 



- ph 



= -AfGSTr [Qt + Qt 



(34) 



where 5 = Ei/u>(Q) is the ratio of the Coulomb exchange 
energy: 



E 1 = — 

TTK 



d 2 q± dq z 

q 2 ± + q 2 z 



\x 2 (q z )\ 



2 p -9i*«/2 



1 



(35) 

and the energy of valley-mixing phonon. Obviously <5 > 
(Ei > 0) and the anisotropy is the easy-plane. The above 
result is valid when Eq <C u{Q)- Otherwise we expect 
this result to be valid qualitatively. The constant (|35| ) in 
the limit I <C Ih reduces t o (H ). 

The sum of polaron ( |29| , |33| ) and the next order E\ 
energies gives the easy-plane anisotropic energy: 



ph 



E ep — G 



J [-Tr(Qr+)Tr(Qr-) 
+ (1 - <5)Tr (Qt+Qi 



2tt ' 



(36) 



We conclude that the JTE is energetically favorable and 
it gives rise to the easy- plane valley anisotropy d36|). 

The easy-plane anisotropy changes dispersion of col- 
lective excitations. Let us consider this effect in the limit 
of vanishing valley symmetry breaking t and Zeeman h 
anisotropics. We add the easy-plane term ( |36|) into the 
Landau-Lifshitz equation (pq): 



Q d t Q 



~Tr(Qi 



+ ) 
(37) 

This equation can be linearized in the vicinity of Q = 
N and describes two 'spin'-modes (unchanged from the 
symmetric case) and one acoustic mode with dispersion 
lu = \/G5E~iq in analogy to the bilayer caseil. 



Summing it up with ( |2^ ) we get the total electron-phonon 
energy in the second order of electron-phonon interac- 
tion. Here we neglect small terms of the order (Q^)~ 4 _due 
to derivative jump of wave function xi z ) & t interfaced At 
v = 1 the energy (|2^^) is isotropic in spin- valley space 
as the two expressions are canceled exactly. If we neglect 
the Coulomb interaction corrections to electron-phonon 
interaction we find that the degeneracy in isospin direc- 
tion is not lifted. Physically, polaron energy is a single 



Phase diagram 

The anisotropic bivalley energy in the uniform state is 
i diagonal matrix element of the anisotropic Hamiltonian 
||) expressed in terms of the order parameter Q: 

E an /N = -t Tt(Qt x ) - h Tv(Qa x ) + 
G [(1 - 5)Tt {Qt + Qt~\ - Tr(Qr + )Tr(Qr-)l , (38) 
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Phase diagram in the v 


= 25 = 


case. F,S,C 



magnetic, spin-singlet and canted antiferromagnetic phases. 

where h = \g e \/iBH. The order parameter Q can be 
parameterized by six {y = 1,3) or eight (y — 2) angles. 
Diagonal matrix elements are real despite the fact that in 
an external magnetic field there is no time reversal sym- 
metry. Therefore the ground state can be chosen real, 
produced from the reference state by SO (4) sub-group 
rotations. This sub-group has 6 parameters with two of 
them falling into the denominator sub-group (y = 2). 
One of the remaining four angles corresponds to global 
rotation of all spins and is fixed by the magnetic field 
direction z. Thus, the ground state differs from the ref- 
erence state by just three rotations {v = 2). 

First we consider the case v = 1. The ground state 
is 17#(T +1), where U = P{d)R(0). Matrix R rotates 
the |-spin component in valley space and then matrix P 
rotates the ±- valley components of the resulting state in 
spin space by angle ±??. Substituting Q(-d } 9) into j3^ ) 
we obtain the energy per one electron: 



E^ n = — h sin d — t sin ( 



GS cos 2 1 



(39) 



The minimum of this energy is reached at i? = n/2 and 
6 = Tt/2 corresponding to the ferromagnet phase in spin 
and valley spaces. The case v = 3 is identical to v = 1 
due to electron-hole symmetry. 

In the case v = 2 the ground state is U^{\ +1,1 —1) 
where U = P(<S)R(9-[,0{). Matrix R rotates the f, i-spin 
components by angles 6f t i in the valley space and then 
matrix P rotates the ±l-valley components by angles ±i? 
in spin space. We find Of = Substituting Q{&, 9) into 
(ph we find for 5 = 0: 



E ln 

T 



-2t cosi? sin 9 — 2/isini9cos# 



+G(1 - 5) I ~ cos 2 9 - cos 2 i?J - G5 cos 2 d sin 2 0. (40) 

The last term represent easy-plane anisotropy and anti- 
ferromagnetic exchange interaction between spins in two 
valleys. The phase diagram in the case 8 = is shown 
on FigJ|. Here F is the spin ferromagnetic and the valley 
singlet phase ($ = w/2, 8 = 0), S is the spin singlet and 
the valley ferromagnetic phase (i? = 0, = 7r/2), C is the 
canted antiferromagnetic phase in spin and valley spaces. 



At S = we find two lines of continuous phase transi- 
tion between F and C phases: (h— G)(h— G/2) = t 2 , and 
between S and C phases: (t + G){t + G/2) = h 2 . Small 
positive 5 transforms the S-C phase transition into a dis- 
continuous first order one. C phase ends at some critical 
h c (S) and the direct S-F first order transition occurs at 
h > h c (5). For 8 > 0.33 C phase disappears and there is a 
line of direct first order S-F transition: t = h—G(3—S) /4. 

At typical experimental magnetic field B ~ 5 — 10T 
Zeeman energy in— Si is few times greater then valley 
splitting t - 2ATii To our knowledge the deforma- 
tion potential at large phonon momentum is unknown. 
Therefore we could only guess the electron-phonon en- 
ergy G ~ 0.2 — 2A. Probably to observe phases of Fig.|| 
in Si one would need to artificially lower the effective 
electron ^-factor. 



Anisotropic energy of skyrmion 

Many experiments found activation gap for diagonal 
conductivity in QHE systems. For symmetry broken 'fer- 
romagnet' systems theory predicts that conductivity is 
mediated by charged topological textures - skyrmions, 
with activation energy being determined by large ex- 
change constant E\ iRjj)- Contrary to this experiments 
find much lower gapsET One possible explanation is that 
due to intrinsic inhomogeneity - defects or long range 
slow variation of electrostatic potential - skyrmions of 
opposite topological charges are already present in the 
system and experiments measure either depinning acti- 
vation energy or skyrmion mobility activation energy. In 
both cases gap is determined by the small anisotropic 
part of skyrmion energy. In this section we find this 
anisotropic energy for a bivalley skyrmion. 

Belavin and Polyakov found jthe skyrmion solution in 
the case of S 2 order parameter!! Their solution is read- 
ily generalized for general Grasmanian order parameter. 
The non-homogeneous order parameter that represents 
one skyrmion in symmetric case with \Q\ = 1 is given by: 



Qbp(z) = N 



1 



R 2 + \z\ 



-R 2 \vf)(v f \ zR\v f ){v e 
zR\v e )(vf\ R 2 \v e )(v e 



(41) 

where z = x + iy, R is the radius of skyrmion core and 
\vf) and \v e ) are the two vectors of size v and 4 — v. 
We choose |w/ e ) = (1,0. .0) in (fi"l|). The skyrmion order 
parameter has to be rotated by a homogeneous matrix 
U calculated in the previous section in such a way that 
the order parameter far away from the skyrmion center 
minimizes the anisotropy energy. In addition to this ro- 
tation we have to allow the global rotations W from the 
denominator sub-group S that transform the skyrmion 
order parameter @: Q(r) = UWQ B p(z)W+U+ . For 
the case v — 2 the matrix W can be parameterized by 
seven angles: 



W 



_(W f 

W e 



(42) 
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where 



W ef = 



COS 



— sin 



sm%-e^- Q ^ 



cos%-e J ^^- c 



(43) 

The additional seventh parameter angle of the denomi- 
nator sub-group rotates the coordinates: z — > e* 77 z. We 
find explicitly that the skyrmion anisotropic energy does 
not depend on phases j e , 7/ and 77 whereas a e = and 
df = 7r correspond to the energy minimum. Important 
are two angles /3 e and (3f which rotate the order param- 
eter inside the core of skyrmion. Besides these angles 
the skyrmion energy depends on conformal parameters 
of the BP solution. We consider here only the simplest 
case of skyrmion topological index: \Q\ = 1 ( [fl| ) where 
only one such parameter is R. There are different spa- 
tial integrals that we encounter in the calculation. One 
integral is logarithmically large: 









(44) 



and we neglect all others. In this way we find Zeeman 
energy of skyrmion: Eg — K(R)Z, where 

Z = 2h (2 sin d cos 9 - cos i?(sin/3/ + sin (3 e ) sin 0) (45) 

The bare valley splitting energy of skyrmion reads: 
E s t k = K(R)T, where 

T = 2t (2cosi?sin6> - sini?(sin/3/ + sin/3 e ) cos 9) (46) 

And finally the Jahn- Teller energy of skyrmion reads: 
G s k — K{R)Q, where 



G = G 



- \ + 3 cos 2 § + i cos 2 9 cos(/3 e + (3 f ) 



- — cos 
2 



2 cos 2 <d 

2 



cos(/3 e - f3 f ) . (47) 



The total energy of a skyrmion includes also a direct 
Coulomb energy due to repulsion of an additional charge 
distribution n(r) inside the skyrmion core (p2|): 



2n\r — r' 



154 ^=^f ^ 



The minimum of the total anisotropic skyrmion energy 
(H,|||,|Z|): £ sk = Z + T + G, over the two parameters e 
and f3f was found numerically and is denoted as: £^ in - 
Next, we find a minimum of the total skyrmion energy 
including @: E sk = K(R)£^ ln +£ C /R, with respect 
to the skyrmion radius R: 



A 



E X \Q\-E Q 



- ( f sk f 2 lotr-^i- 

2 I ^min^C 10 § csk 



1/3 



(49) 



This formula is valid in the limit £^ k in -C Ei . The result- 
ing anisotropic part of a skyrmion gap is shown on the 



FIG. 4: Anisotropic part of skyrmion energy gap at v = 2 
and S — 



Fig.||. Note that the prominent minimum of the skyrmion 
anisotropy energy gap coincide with the phase transition 
line C-F from the phase diagram Fig|3[ 

This skwrmion anisotropic energy is similar to the bi- 
layer caseE3 with t being the hopping constant between, 
the two layers. In the experiment on GaAs bilayerlia 
a profound reduction of the thermal activation gap has 
been found in some interval on the v = 2 line. 

In the case v = I we parameterize general rotations 
from the denominator sub-group by four angles: 



cos — , sin — cos are , sin — sin ae ' 
2' 2 '2 



(50) 



The skyrmion energy does not depend on the phases Ai^ 
and a = tt/2 in the energy minimum. We find Zeeman 
energy Z = h(l — cos f3), the valley splitting energy T = 
2tcos 2 (/3/2) and the easy-plane Jahn- Teller energy g — 
GS (I + cos (3). The minimum of these skyrmion energies 
over the parameters (3 and a: £^ in = 2 min (t + G5, h), 
determines the charge activation gap A in Eq.([l9|). The 
case v = 3 is identical to the case v = 1. 



Conclusion 

We have shown that the valley degeneracy for Si(IOO) 
MOSFET can lead to Jahn Teller effect in strong mag- 
netic field with a longitudinal lattice deformation inside 
the width of proper quantum well. This deformation has 
an atomic periodicity corresponding to the momentum 
difference between the two valley minima in the direction 



8 



perpendicular to the plain. Integer fillings v = 1,3 and 
v = 2 are treated in different ways. For v = 1, 3 the Jahn 
Teller effect alone can not remove the valley degeneracy 
and Coulomb interaction is crucial. The level splitting 
here is small if the Coulomb energy is small compared to 
the Debay energy. For v — 2 case Jahn- Teller effect re- 
moves the valley degeneracy via an easy-plane anisotropy 
similar to that found in the 2D bilayer system in strong 
magnetic field. Phase diagram depending on different 
physical parameters was established and anisotropic en- 
ergy of topological skyrmion like texture was calculated. 

This work was supported by INTAS grant 99-01146 
and RFBR grant 01-02-17520a. Author express their 
gratitude to V.M.Pudalov and V.T.Dolgopolov for use- 
full discussions. 



(G 00 (z'Or> + (£)Goi(^) + Gio^'OMOGoofc)) , (52) 

Expanding fi(£) around the center point of the diagram: 
z = (z + z')/2: n±(0 = fl+(z )±(£-z ) e M „V^(z ), 
we evaluate the energy 6E (^2|) as (z = r): 

J d 2 z (l - Jf£) y(M)e-H 2 /2 / e^^Hf) ^ 

(53) 

The front integral is evaluated for the Coulomb interac- 
tion and equals: (Ei + Eq)/2. 



Appendix A 



Appendix B 



Here we calculate the_second diagram on Fig.[l] which 
was overlooked in RefEfl. It represents the change of 
Hartree-Fock exchange energy in a presence of spin tex- 
ture. The spatial part of the electron Green function in 
Landau gauge is: 

Gos ( Zj Z >) = l Z ~Zy e -\ Z - Z '\Z/4+i( X+X ')(y-y>)/2^ (51) 

V2 s s! 

where z = x + iy, s is the Landau level and G s o = Gq s . 
The bottom part of Figjljb features Goo Green function 
whereas the upper part has Gio and Goo Green functions. 
Taking the two frequency integrals we find: 



SE = —= 



1 [d 2 zd 2 z'd 2 Z 



2tt 2tt 2tt 



G 00 (z,z')V( Z -z') 



Eq. (^5|) is prooved using a principal bundle of complex 
Stiffel manifold V£ k = SU{n)/SU{n - k) - that is de- 
fined as a manifold of k orthogonal complex vectors in 
n dimensional complex linear space - over Grassmanian 
manifold G^ k with a layer U(fc). The exact map sequence 
for this bundle is 



Stiffel manifold has a property nj(V r 



c\ - 



*i{Vg)- (54) 
at j < 2(n - 



kfB : proven using the principal bundle SU(n) over V^ k 
with a layer SU(n — k). Eq.(54|) means that 7r 2 (G^ fe ) = 
7ri(t/(fc)). The Lee group U(fc) is a product of SU(fc) 
and U(l) groups with the fundamental homotopy group: 

7ri(c7(fc)) = ■K X {su{k)) + m(u(i)) = zm 
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